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ABSTRACT 

Aims. To constrain the mass-to-radius ratio of isolated neutron stars, spin-phase resolved X-ray spectroscopic analysis is performed. 
Methods. The data from all observations of RBS 1223 (IRXS J130848.6+212708) conducted by XMM-Newton EPIC pn with the 
same instrumental setup in 2003-2007 were combined to form spin-phase resolved spectra. A number of complex models of neutron 
stars with strongly magnetized (Spoie ~ 10'^ - lO''* G) surface, with temperature and magnetic field distributions around magnetic 
poles, and partially ionized hydrogen thin atmosphere above it have been implemented into the X-ray spectral fitting package XSPEC 
for simultaneous fitting of phase-resolved spectra. A Markov-Chain-Monte-Carlo (MCMC) approach is also applied to verify results 
of fitting and estimating in multi parameter models. 

Results. The spectra in different rotational phase intervals and light curves in different energy bands with high S/N ratio show a high 
complexity. The spectra can be parameterized with a Gaussian absorption line superimposed on a blackbody spectrum, while the light 
curves with double-humped shape show strong dependence of pulsed fraction upon the energy band (13% -42%), which indicates 
that radiation emerges from at least two emitting areas. 

Conclusions. A model with condensed iron surface and partially ionized hydrogen thin atmosphere above it allows us to fit simul- 
taneously the observed general spectral shape and the broad absorption feature observed at 0.3 keV in different spin phases of RBS 
1223. It allowed to constrain some physical properties of X-ray emitting areas, i.e. the temperatures (Tpi ~ 105 eV, Tp2 ~ 99 eV), 
magnetic field strengths (Spi a Sp2 ~ 8.6 x 10'^ G) at the poles, and their distributions parameters (ai ~ 0.61, 02 ~ 0.29, indicating 
an absence of strong toroidal magnetic field component). In addition, it puts some constraints on the geometry of the emerging X-ray 
emission and gravitational redshift (z = 0.16!jj j;^) of RBS 1223. 
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1. Introduction Much more information may be extracted by studying spin- 

^""^ phase resolved spectra with high signal to noise ratio and fitting 

^ , The study of thermally emitting, radio-quiet, nearby, isolated ^^em with model spectra of radiation emitted from highly mag- 

L| ' neutron stars (INSs) may have an important impact on our un- netized INS surface layers 
' derstanding of the physics of neutron stars. Observations and 
modeling of thermal emission from INSs can provide not only 
5^ information on the physical properties such as the magnetic field, RBS 1223 (IRXS J130848. 6+212708) was originally dis- 

temperature, and chemical composition of the regions where this covered as a soft X-ray source during the ROSAT All-Sky 

radiation is produced, but also information on the properties of Survey by ISchwope et"al] |l999). It shared common charac- 

matter at higher densities deeper inside the star. teristics with other members of the small group (so far 7 

In particular, measuring the gravitational redshift of an iden- discovered by ROSAT) of thermally emitting and radio-quiet, 

tified spectral feature in the spectrum of thermal radiation emit- nearby INSs, traditionally dubbed XDINS (fror n "X-ray dirn 

ted from the INS surface/atmosphere may provide a useful con- INS") or Magnificent Seven (see, e.g., reviews bv Haber fcoOTl : 

straint on theoretical models of equations of state for super- Mereghetti 2008; Turolla 2009, and references therein): soft 

dense matter. Indepe ndent of the estimate of the INS radius (e.g. spectra, well described by blackbody radiation with tempera- 

iTriimper et al.ll2004) from the thermal spectrum of an INS, it tures 60- 120 eV, no other spectral features, no association with 

may allow us to directly estimate the mass-to-radius ratio. SNR, no radio emission, no X-ray pulsations, very high X-ray to 

optical flux ratio. A nature of these sources as old INS reheated 

* Based on observations obtained with XMM-Newton, an ESA sci- by accretion from the interstellar medium or young cooling stai's 

ence mission with instruments and contributions directly funded by seemed possible. Meanwhile, intensive X-ray observations with 

ESA Member States and NASA mainly XMM-Newton and also with Chandra, as well as opti- 
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cal/UV observations of most probable counterparts, have revised 
this picture in parts and has provided intriguing physical insight. 

X-ray pulsations have been found in six of these objects, 
with periods clustered at 3 - 11 sec, and period derivatives have 
been measured for five of them. In a classic P - P diagram 
the XDINSs are found intermedi ate between radio pulsars and 
magnetars (e.g.. Mere ghettill2008h . The inferred magnetic field 
strengths are above 10'^ G. 

Imaging CCD-spectroscopy with XMM-Newton has uncov- 
ered absorption features in at least three, likely in six stars. At 
current energy resolution (FWHM ~ 50-150 eV between 0.2- 
2.0 keV) they can be, formally, well fitted as a Gaussian ab- 
sorption lines, which are usually connected with ion cyclotron 
lines. Their interpretation is not unique, magnetically shifted 
atomic transitions or electron cyclotron resonances are debated. 
For some of the mo dels, the inferr ed magnetic field strengths are 
again above 10'^ G (lHaberlll2007h . 

Among them, RBS 1223 is a special case, some sort of out- 
lier, qualitatively different from the typical pattern observed in 
the sample of XDINSs. Namely, the rotational phase folded light 
curve has a double-humped shape, with largest pulsed fractiorQ 
(~ 19% in 0.2-1.2 keV energy range, see further). Moreover, the 
separation of maxima by less than 180 degrees, significantly dif- 
ferent count-rates of minima, and the variation of the blackbody 
apparent temperatures between 80 and 90 eV over the spin cycle 
are evidence for, at least, two emitting area s with some tempera- 
ture d istribution over neutron star surface (ISchwope et al.ll2005l 
12001) . 

It is worth to note that the abovementioned observed absorp- 
tion feature at ~ 0.3 keV in the sp ectrum of RBS 1223 has the 
highest equivalent width (~ 0.2keV. [Schwope et al.l2007l) among 
all XDINSs. Moreover, spectral analysis of the avera ge XMM- 
Newton spectrum of RBS 1223 ( Schwop e et alj|2007 h based on 
the two first observations (see Table[T]) showed the possible pres- 
ence of a second feature in the X-ray spectrum. Its existence, 
however, is not uniquely proven because of some remaining cal- 
ibration uncertainties and the not well-defined continuum at high 
energies because of the lack of photons (lower S/N ratio). 

Meanwhile, new observational sets (see Table [1]) are pub- 
licly available and a number of new detailed models of emergent 
spectra of highly magnetized INSs are d eveloped and available 
dHo et al.ll2009tlSuleimanov et al.ll2010ah . 

Preliminary analysis of the phase averaged spectrum of RBS 
1223 was performed by Perez-Azorm et al. (2006b). They men- 
tioned that there is a possibility of good fits with quadrupo- 
lar magnetic fields, although not excluding a condensed surface 
model with h ydrogen atmosphere , including vacuum polariza- 
tion effects (van Adelsberg & Lai 2006). 

To explain unusual observed properties of RBS 1223 
Suleimanov et al. (2010a, hereafter Paper I) studied various lo- 

' The pulsed fraction, used in this paper, is defined as: 
PF = , 

where CR is the count rate. Note that this definition of a pulsed frac- 
tion may be misleading in the case of complex shaped light curves, a 
peak/minimum flux may appear at different phases in different energy 
ranges. Instead, the following quantity, i.e semi-amplitude of modula- 
tion, might be an appropriate descriptor of the pulsed emission: 

_ I], ICR, - (CR,>| 
2,CR, 

where CR, is a count rate per phase bin. 
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Fig. 1. Observed count rates of RBS1223 in the energy band 0.16 
- 2.0 keV in different XMM-Newton EPIC pn observations. 

cal models of the emitting surface of this INS. They considered 
three types of models: naked condensed surfaces, semi-infinite 
partially ionized hydrogen atmospheres with vacuum polariza- 
ti on and partial mode con version taken into account (see details 
in lSulermanov et al.l2009l) and such thin atmospheres above con- 
densed iron surface. They also created a code for modeling inte- 
gral phase resolved spectra and light curves of rotating neutron 
stars. This code takes into account general relativistic effects and 
allows one to consider various temperature and magnetic field 
distributions. Analytical approximations of the considered three 
types of local spectra were used for such modeling. It was qual- 
itatively shown that only a thin model atmosphere above a con- 
densed iron surface can explain the observed equivalent width of 
the absorption feature and the pulsed fraction. 

In this paper we use the code developed in Paper I to perform 
a comprehensive study of co-added high-S/N and phase-resolved 
XMM-Newton EPIC pn spectrum of RBS 1223, despite a possi- 
ble small variations of brightness shown by this INS (see Fig.[Tli. 

2. Observations and data reduction 

RBS 1223 has been observed many times by XMM-Newton 
(Table [Hi. Here we foc us on the data collected with EPIC pn 
(den Her der et al.ll200ll) from the 12 publicly available observa- 
tions, with same instrumental setup (Full Frame mode, Thinl 
filter, positioned on axis), in total presenting about 175 ks of ef- 
fective exposure time. 

The data were reduced using standard threads from the 
XMM-Newton data analysis package SAS version 10.0.0. We re- 
processed all publicly available data (see Table [1]) with the stan- 
dard metatask epchain. To determine good time intervals free of 
background flares, we applied filtering expression on the back- 
ground light curves, performing visual inspection. This reduced 
the total exposure time by ~30%. Solar barycenter corrected 
source and background photon events files and spectra were 
produced from the cleaned SINGLI^ events, using an extrac- 
tion radius of 30" in all pointed observations. We extracted also 
light curves of RBS 1223 and corresponding backgrounds from 
nearby, source free regions. We then used the SAS task epiclc- 
corr to correct observed count rates for various sorts of detector 
inefficiencies (vignetting, bad pixels, dead time, effective areas, 

- See XMM-Newton Users Handbook. 
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Fig. 2. Energy-Phase image of RBS 1223 combined from dif- 
ferent observations. Rotational phase-folded light curve in the 
broad energy, 0.2-2.0 keV, band (left panel) and the phase aver- 
aged spectrum (bottom panel) are shown. 

Table 1. XMM-Newton EPIC pn observations of RBS 1223 



Obs. ID 


Obs. Date 


Exposure 


Effective exposure 




MID 


ksec 


ksec 


0157360101 


52640.4271023 


28.910 


25.966 


0163560101 


53003.4674787 


32.114 


26.436 


0305900201 


53546.0889925 


16.806 


11.443 


0305900301 


53548.0919318 


14.812 


12.682 


0305900401 


53566.4161765 


14.814 


10.574 


0305900601 


53745.8803029 


16.845 


14.702 


0402850301 


53894.9783016 


7.419 


4.761 


0402850401 


53902.9489095 


8.421 


5.558 


0402850501 


53913.2292603 


12.517 


2.801 


0402850701 


53921.2857528 


10.411 


6.366 


0402850901 


54096.6780925 


9.315 


8.011 


0402851001 


54262.6501350 


10.920 


8.462 



etc.) in different energy bands for each pointed observation (see 
Fig.D. 

For each registered photon the corresponding rotational 
phase was computed according to the phase coherent timing 
solution provided by Kaplan & van Kerkwiik (2005). We pro- 
duced spin phase-folded light curves in different energy bands 
and spectra corresponding to different phase intervals for each 
pointed observation. Finally, the latter ones were co-added for 
the same intervals of phases to create the combined phase- 
resolved spectra of RBS 1223 (Fig.|2]i. 

It should be noted that spectral responses and effective areas 
for those 12 different observations were almost undistinguish- 
able. 



3. Data analysis 

First, we have fitted the phase-averaged, high signal-to-noise ra- 
tio spectral data, collected from different observations (Fig. [3] 
and TablelJI, using a combination of an absorbed blackbody and 
a Gaussian absorption line multiplicative component model (in 
XSPEC phabs*bbodyrad*gabs). 



In the spectral energy range of 0.16-2.0 keV we obtained a 
statistically acceptable fit with the reduced = 1.1. 

We performed also a fit with the abovementioned spectral 
model for four phase intervals, including maxima and minima 
of the double-humped phase folded light curve, simultaneously. 
The fitted parameters are presented in Table|2] It shows a depen- 
dence of the blackbody apparent temperature and the Gaussian 
absorption feature center upon rotational phase. Next, we di- 
vided the broad energy band into four energy regions (0.16- 
0.5, 0.5-0.6, 0.6-0.7, and 0.7-2.0 keV) and constructed rota- 
tional phase folded light curves (Fig. |4]l. It is also noteworthy 
that the higher the considered spectral energy range, the larger 
the pulsed fraction. These spin phase-folded light curves may 
serve for rough estimates of the viewing geometry and physi- 
cal characteristics of emitting areas of RBS 1223 (see below). In 
particular, the obtained parameters of the fits at different phases 
(Table |2]l support a model with hot areas around the magnetic 
poles, because the energies of line centers and the blackbody 
temperatures are larger at the peaks. 

In view of this result, we implemented into the X-ray spec- 
tral fitting package XSPEC a number of new highly magnetized 
INS surface/atmosphere models developed in Paper I. They are 
based on various local models and compute rotational phase de- 
pendent integral emergent spectra of INS, using analytical ap- 
proximations. The basic model includes temperature/magnetic 
field distributions over INS surface, viewing geometry and grav- 
itational redshift. Three local radiating surface models are also 
considered, namely, a naked condensed iron surface and partially 
ionized hydrogen model atmospheres, semi-infinite or finite on 
top of the condensed surface. Here we have reproduced the es- 
sential part of the basic model (for details and further references, 
see Paper I). 

To compute an integral spectrum, the model uses an analyt- 
ical expression for the local spectra: a diluted blackbody spec- 
trum for both semi-infinite and thin models of a magnetized at- 
mosphere with one absorption feature; 



hia) = D Be(T) <p(a) exp(-T), 



(1) 



where a is the angle between radiation propagation direction and 
the surface normal, and (pia) represents the considered angular 
distributions of the specific intensities for different cases of at- 
mosphere models. We have represented ip{a) by the following 
three models: 



<P(a) 



1, 

0.4215 -H 0.86775 cosa. 

1 - 0.2 cos- Q-c 

if 



cos a > cos fff 



1 - cos2 ' " ( and E, 
0.2, if cos a < cos and £, 
1 for other energies. 



c, < £ < 4Ec, 



<E < 4Ec, 



Here, the cases a, b, and c correspond to the models of isotropic 
emission, electron scattering directivity pattern, and a finite hy- 
drogen atmosphere layer above a condensed iron surface, respec- 
tively; cosffc = |^|(£'c, /£c - 1)] ^ is a parameter of angular 
distribution for the case of a thin atmosphere in the energy range 
Ea < E < 4Ec', ^ci = fi ZeB/niiC and Ec.e - lieB/mgC are the 
ion and electron cyclotron energies, Zip e - fi ^jA^Te^n^|m^ is the 
electron plasma energy, Z and nii are the ion charge number and 
ion mass, is the electron number density, B the magnetic field 
strength. 



(2) 
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Table 2. Fitting results by model phabs*gabs*bbodyrad (r = 0.612 + 0.08, cr = 0.155 + 0.008, A^h = (1-8 ± 0.0035) x lO'^cm^^) 



Phase interval 



kT 

(keV) 



Line center 
(keV) 



Remarks 



0.20-0.35 
0.35-0.65 
0.65-0.80 
0.0-0.2,0.8-1.0 
0.00-1.0 



0.084±0.001 
0.088±0.001 
0.083±0.001 
0.088±0.001 
0.086±0.001 



0.25+0.01 
0.26±0.01 
0.24±0.01 
0.29±0.01 
0.27+0.01 



First minimum 
Secondary peak 
Second minimum 
Primary peak 
Phase averaged 




Fig. 3. Phase-averaged X-ray spectrum (left panel) and phase-folded light curve (right panel) in the broad energy band of 0.16-2.0 
keV of RBS 1223 combined from 12 pointed XMM-Newton EPIC pn observations. 




Fig. 4. X-ray spectra including primary and secondary peaks, first and second minima, and phase-folded light curves in different 
energy bands of RBS 1223 combined from 12 pointed ZMM-A^evvfon EPIC pn observations. Fitted absorbed blackbody with gaussian 
absorption line models {phabs'*bbodyrad*gabs XSPEC) to the spectra and two sinusoids to the phase-folded light curves are also 
shown (for details, see text and Table|2]). 

Dotted vertical lines are indicating phase intervals used for extraction and fitting of spectra shown in the left panel. 



and 

£c, =£c (1+3(1 -cos a)^/2)^ (3) 

The optical depths t - r° exp^- ^"^'^;"-* j and the widths 

o"iine of the absorption features are considered identical for all 
local spectra, but the center of the line depends on the local mag- 
netic field strength. In the case of a semi-infinite atmosphere this 
absorption line represents the blend of the proton cyclotron line 
and nearby b-b atomic transitions in neutral hydrogen. The line 



center corresponds to the proton cyclotron energy E^ ^. In the 
case a thin atmosphere the absorption feature was represented 
by a half of a Gaussian line. It means that t" = at £ < E\\ag. 
This half of the Gaussian line represents a the complex absorp- 
tion feature, which includes a broad absorption feature from the 
emitting condensed iron surface transmitted through a thin hy- 
drogen atmosphere and the proton cyclotron line plus b-b atomic 
hydrogen transitions (see details in Paper I). The center of this 
line corresponds to the ion cyclotron Ec,i for completely ionized 
iron. In both cases, the absorption details are better represented 
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by sums of two or even three Gaussian lines. But we have chosen 
to use only one (or half) Gaussian line to reduce the number of 
fitting parameters. The local spectra of the condensed iron sur- 
face is approximated by simple step functions (see Paper I). The 
dilution factor D cannot be found from the fitting independently, 
and we take D - 1 for all local spectra. 

The blackbody temperature and the magnetic field distribu- 
tions in the two emitting areas around the magnetic poles were 
presented analytically bv lPerez-Azorm et al. I (l2006ah : 



cos 



- pl,2 



cos^ u + a[ 2 sm 



and 



B - Bpi_2 ^Jcos^ 6 + fli_2 sin 6, 



(4) 



(5) 



where T-p\ 2 and B^i 2 are the temperatures and magnetic field 
strengths at the poles and a\ 2 their distribution parameters and 
Jniin is the minimum temperat ure reached on the surface of the 
star (we chose « O.STp, like in lPerez-Azorm et"a n i2006ah . Here 
9 is the magnetic colatitude. 

The parameters 012 are approximately equal to the squared 
ratio of the magnetic field strength at the equator to the field 
strength at the pole, a\ 2 ~ (Beq/Bpi,2)^- Using these parameters 
we can describe various temperature distributions, from strongly 
peaked (a » 1) to the classical dipolar (a - 1/4) and homoge- 
neous (fl - 0) ones. 

Finally, the total observed flux at distance d from the INS 
is summed over the visible surface. A local spectrum at photon 
energy E from the unit surface is computed taking into account 
viewing geometry, gravitational redshift, and Ught bending ef- 
fects (e.g., Poutanen & Gierliriski 2003h : 



/e t/^c/y = Cnoim (1 + z) lE'ia) cosa sinydipdy. 



(6) 



where Cnoim = R^/d^ is a normalization constant and y is the 
angle between the radius vector at a given point and the rota- 
tion axis. The observed and emitted photon energies are related 
as £■ = E'{1 + z), where z is gravitational redshift. The an- 
gle a between the emitted photon and normal to the surface in 
the local reference frame depends on z, the local surface posi- 
tion (on azimuthal and colatitude angles (p and y), and the angle 
between rotation axis an d line of sight and is computed us- 
ing the approximation of iBeloborodo vl ()20Q2l). The summation 
is performed separately for both regions around magnetic poles, 
because we allow the poles to be not precisely antipodal (one 
of them can be shifted by angle k relative to the symmetric an- 
tipodal position). Therefore the poles may have different angles 
between rotation and magnetic axes 0b 1,2- 

As shown in Paper I, a model with a thin hydrogen atmo- 
sphere above the condensed iron surface with a smooth temper- 
ature distribution over the neutron star surface may describe very 
well the observed physical properties of RBS 1223. 

This basic model has a number of input parameters depend- 
ing on the inner atmosphere boundary condition of an INS (con- 
densed iron surface or blackbody, temperature and magnetic 
field distributions over surface, viewing geometry and gravita- 
tional redshift) and angular distribution of the emergent radiation 
(isotropic or peaked by a thin atmosphere above a condensed 
surface). 

Combined spectra of RBS 1223 in 20 phase bins (Fig. |2]l 
were considered as different data sets during simultaneous fitting 
with the absorbed abovementioned model. All input parameters 
were free and linked between those data groups. 



Given the large number of the free parameters in the model 
we performed a preliminary analysis in order to get rough esti- 
mates of (or constraints on) some of them. From the observed 
double-humped light curve shape in different energy bands (see 
Fig. nil, it has been akeady clear that two emitting areas have 
different spectral and geometrical characteristics (e.g., the rela- 
tively cooler one has a larger size). Moreover, from the peak-to- 
peak separation in this double-humped light curve we locate the 
cooler one at an offset angle of k with respect to the magnetic 
axis and azimuth (ISchwope et alJl2005h . 

First we performed formal spectral fitting with the simplest 
model, i.e. absorbed blackbody with multiplicative gaussian ab- 
sorption line including phases of maxima (see Table |2]i in the 
Ught curve (Figs. [3] and |4|. It is clear, that at these phases X- 
ray emission is mostly dominated by emitting areas at the mag- 
netic poles, where temperatures also have maxima (see formula 
of temperatur e and magnetic field dependence upon polar an- 
gle, Perez-Azorfaeral]2006i and Paper I). Secondly, some con- 
straints on magnetic field strengths at the poles may be used on 
the base of period and its derivative values assuming magnetic 
dipole breaking as a main mechanism of the spin down of RBS 
1223. 

Having these two general constraints on temperatures and 
magnetic field strengths at the poles we simulated a large number 
of photon spectra (absorbed blackbody with gaussian absorption 
line and different models considered in the Paper I) folded with 
the response of XMM-Newton EPIC pn camera, taking into ac- 
count also the interstellar absorption (for parameters see Table|2] 
and using a characteristic magnetic field strength of B~ 3.4x10'^ 
Gauss). 

The predicted phase-folded light curves in four spectral 
ranges: 0.16-0.5 keV, 0.5-0.6 keV, 0.6-0.7 keV, 0.7-2.0 keV and 
with free parameters of viewing geometry and gravitational red- 
shift, normalized to the maximum of the brightness, we cross- 
correlated with observed ones (Fig.|4]i. We infer some constraints 
on the parameters from the unimodal distributions of them when 
cross-correlation coefficients were exceeding 0.9 in the men- 
tioned four energy bands simultaneously and used them as initial 
input value and as lower and upper bounds for fitting purposes. 
For example, gravitational redshift cannot exceed 0.3 (it is not 
possible to obtain the observed PF at larger z due to strong light 
bending) or the antipodal shift angle must be less than 25° (due 
to the observed phase separation between the two peaks in the 
light curve). It is also evident that the sum of the inclination an- 
gle of the line of sight and magnetic poles relative to the rota- 
tion al axis are already constrain ed by the light curve class (see, 
Pou tanen & Beloborodovll2006l class III) and have the maxi- 
mum effect (e.g., provide the maximum PF) when both are equal 
to 90°. It should be noted that these angles do not have a large 
influence on the fitting, and the only important issue is a range 
of values when two maxima are observed. 

Having abovementioned crude constraints and input values 
of free parameters we performed fitting with the models imple- 
mented in XSPEC package of combined spectra of RBS 1223 in 
20 phase bins (Fig.|2]i simultaneously, i.e. each of those phase re- 
solved spectrum considered as different data sets with the linked 
parameters to the others, and the only differences were phase 
ranges, which were fixed for an individual phase resolved spec- 
trum. 

The fitting was successful, with C-statistic value 2937 with 
2159 degrees of freedom. These parameters are presented in 
Tabled 

In order to assess a degree of uniqueness and to estimate 
confidence intervals of the determined parameters, we have ad- 
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Table 3. Simultaneous fitting results of combined, phase resolved X-ray spectra of RBS 1223 with different spectral models 



Fitted 


Spectral Model 




Parameter* 


iron condensed surface 


BlacKDody 




partially ionized H atmosphere 


electron scattering 


^pi [eV] 


io5.o!y;j 


109 ±4.0 


T^2 [eV] 


99.0 ± 3.0 


106 ± 3.0 


Spi X 10"* [G] 


0.86 ± 0.02 


0.65 ± 0.03 


Sp2 X 10'-* [G] 


0.86 ± 0.02 


0.58 ± 0.02 




0.61 ± 0.11 


0.25 




0.29 ± 0.05 


0.25 


y 


2 76+""' 

^- ' "-0.02 


1.90 ±0.06 


a- [eV] 


225.8+f° 


168.0 ±6.0 


z 


0.15 ±0.02 


0.17 ±0.03 


41 


4.2 ± 0.6 


4.0 ± 0.4 


in 


48.9 ±0.5 


45.6 ±0.5 


en 


90.0 ± 0.5 


90.0 ± 0.5 



Notes. See section[3]for definition of the model parameters 




Fig. 5. ProbabiUty density distributions of parameters (tempera- 
tures, magnetic field strengths at the poles, antipodal shift angle 
and gravitational redshift) by MCMC fitting with the model of 
a neutron star with condensed surface and partially ionized hy- 
drogen layer above it. The most probable parameter value is in- 
dicated by the solid vertical line. Dashed vertical lines indicate 
the highest probability interval (68%, for details see text). 



ditionally performed Markov Chain Monte Carlo (MCMC) fit- 
ting as implemented in XSPEC. In Fig. |5]we present probability 
density distributions of some of them. Note, independent initial 
input values of parameters of MCMC approach converged to the 
same values, in 6 difi'erent chains. 

4. Discussion 

The combined phase-resolved spectra of RBS 1223 can be si- 
multaneously fitted by emergent radiation of a spectral model 
of an iron condensed surface with a partially ionized hydrogen 
atmosphere above it. Formally they can be fitted also by a black- 
body spectrum with proton-cyclotron absorption gaussian line 
and a peaked (typical for an electron scattering atmosphere) an- 
gular distribution of the emergent radiation. In both cases two 
emitting areas with slightly different characteristics are required 
(see TableO. 

It is worth to note, that the resulting fit parameters are very 
similar for different spectral models (see Table O, which is also 
confirmed by an MCMC approach with dififerent input parame- 
ters. 



However, we believe that the emission properties due the 
condensed surface model with a partially ionized, optically thin 
hydrogen layer above it, including vacuum polarization effects, 
is more physically motivated. Moreover, semi-infinite atmo- 
spheres have rather fan-beamed emergent radiation (see Paper 
I and references within) and it seems impossible to combine a 
proton cyclotron line with a pencil-beamed emergent radiation. 

We calculated a set of thin highly magnetized partially ion- 
ized hydrogen atmospheres above a condensed iron surface with 
magnetic field strength B = 8 x 10'"^ G, which is close to 
the value estimated from observations. The observed blackbody 
temperature of the spectra is reproduced at effective tempera- 
tures Teif ~ 7 X 10^ K. Examples of the computed emergent 
spectra with parameters B = 8 x lO'"* G, Tefi = 7 x 10^ K and 
various atmosphere column densities Z are shown in Fig.|6] The 
emergent diluted blackbody spectrum with T = 0.1 keV and 
the absorption Gaussian line parameters, presented in Table |3] 
(t° = 2.8,cr,i„e = 0.226 keV, £iine = 0.24 keV) andD = 0.34 are 
also shown. Unfortunately, from these models we cannot evalu- 
ate the actual atmosphere thickness (the surface density), but we 
can obtain the dilution factor D, which is important to coiTect 
the distance estimation. 

Emission spectra based on realistic temperature and mag- 
netic field distributions with strongly magnetized hydrogen at- 
mospheres (or other light elements) are formally still an alterna- 
tivqj, but it is unphysical because of the absorption line added 
by hand. 

A purely proton-cyclotron absorption line scenario can be 
excluded owing to the equivalent width of the observed ab- 
sorption spectral feature in the X-ray spectrum of RBS 1223. 
Magnetized semi-infinite atmospheres predict too low an equiv- 
alent width of the proton cyclotron line in comparison with the 
observed one. 

This result of the fitting (with the condensed surface model 
with partially ionized, optically thin hydrogen atmosphere above 
it, including vacuum polarization effects) suggests a true ra- 
dius of RBS 1223 of 16 ± Ikm for a standard neutron star 
of 1.4 solar mass, considerably larger than the canonical ra- 

^ Our attempt to fit the combined, phase-averaged spectrum of RBS 
1223 by partially ionized, strongly magnetized hydrogen or mid-Z ele- 
ment plasma model {XSPEC «.?may. iMori & Ho 2007; Ho et al. 2008D, 
as well two spots or purely condensed iron surface models, failed. 
Noteworthy, an acceptable fit is obtained by nsmax model with addi- 
tional, multiplicative gaussian absorption line component (model gabs). 



6 



V. Hambaryan et al.: Phase resolved spectroscopy of RBS 1223 




Photon energy [keV] 



Fig. 6. Emergent spectra of the fitted model of a magnetized at- 
mosphere with a condensed iron surface and a partially ionized 
hydrogen layer above it (see Table |3]l. Different lines are cor- 
responding different atmospheric surface densities. For compar- 
ison purposes a blackbody spectrum with gaussian absorption 
line is also shown in gray (see text). 




Fig. 7. Mass-r adius relations for s everal equations of state (thin 
solid curves, iHaensel et al. I l200l . and a strange star (dashed). 
Thin dash-doted: curves of constant apparent radius Roa - 
Rl ^Jl - IGMjRc^- = 9,13 and 17 km (lTrumpeil2005h . 



dius of 10 km; it is only marginally compatible with the range 
from =s 10 km to ~ 14 km, allowed by modern theoretical 
equations of state of superdense matter (Haensel et al. 200^ 
iHebeler et al.P 2010, and references therein), and indicates a very 
stiff equation of state of RBS 1223 (Fig . [71 for similar results, see 
alsojH 
l2006r 



b et al.ll2007HHeinke et al.l2006i:[Suleimanov & PoutanenI 



Suleimanov et al.ll2010bl) " 



With these estimates of the radius and normalization con- 
stant of the fit (Rjd ~ 0.0247 km pc"'), we obtained for 
the assessment of the distance of RBS 1223 VD650^jqPC - 
V£>/0.34 380+1?, (ggg also lSchwope et al.ll2005h . 



5. Conclusions 

The observed phase resolved spectra of the INS RBS 1223 are 
satisfactorily fitted with two slightly different physical and ge- 
ometrical characteristics of emitting areas, by the model of a 
condensed iron surface, with partially ionized, optically thin hy- 
drogen atmosphere above it, including vacuum polarization ef- 
fects, as orthogonal rotator. The fit also suggests the absence of 
a strong toroidal magnetic field component. Moreover, the de- 
termined mass-radius ratio {{M / Mq)/ (R /km) = 0.087 ± 0.004) 
suggests a very stiff equation of state of RBS 1223. 

These results on RBS 1223 are promising, since we could 
find good simultaneous fits to the rotational phase-resolved spec- 
tra of RBS 1223 using analytic approximations for the above- 
mentioned model implementation. 

More work for detailed spectral model computation will 
be certainly worth to do in the near future and application to 
the phase-resolved spectra of other INSs. In particular, for RX 
J1856.5-3754 and RX J0720.4-3125, including high resolu- 
tion spectra observed by XM M-Newton and Chandra with possi- 
bly other absorption features dHambarvan et alj|2009l : iPotekhinI 
20101) . 
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